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ABSTRACT 


Simple  analytical  methods  for  the  determination  of  the  angle  of 
emergence  from  complex  poles  and  zeros  are  developed.   Several  case 
studies  of  root  loci  emerging  from  complex  zeros  are  analyzed  to  deter- 
mine criterion  for  the  loci  starting  from  complex  pole  to  end  on  com- 
plex zero,  as  a  pre-requisite  study  for  "cancel -compensation"  and  self- 
adaptive  compensators. 
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CHAPTER  I 
INTRODUCTION 

Many  closed  loop  transfer  functions  of  control  systems  such  as 
those  representing  structual  resonances  in  missiles,  fuel  sloshing, 
shaft  twisting,  etc.,  have  a  pair  of  complex  conjugate  poles  near  the 
imaginary  axis.   If  gain  requirements  are  not  too  high,  the  effect  of 
the  roots  near  these  poles  is  often  neglected  because  other  roots 
dominate.   However,  in  the  design  of  high  gain  systems,  these  roots 
near  the  imaginary  axis  may  cross  over  into  the  right  half  s-plane  and 
thus  cause  instability  in  the  system.   Figure  1  shows  the  effect  of 
raising  gain  on  the  open-loop  Bode  diagram. 

Cancel  compensation  would  be  an  effective  solution  to  removing 
these  complex  poles  near  the  jw-axis  if  it  were  not  for  (1)  the  diffi- 
culty in  the  design  accuracy  of  these  compensators,  and  (2)  the  poles 
themselves  may  not  be  fixed,  and  may  move  about  in  an  area  in  the  s- 
plane. 

The  mere  placing  of  complex  zeros  near  the  complex  poles  is  no 
assurance  that  the  root  loci  is  connected  between  these  poles  and 
zeros.   (See  Fig.  2)   When  this  loci  does  close  on  the  adjacent  com- 
plex zero,  the  root  on  this  locus  has  such  a  small  residue  that  it  can 
be  neglected  when  compared  with  the  other  roots. 

The  question  remains:   What  are  the  requirements  for  the  location 
of  complex  zeros  adjacent  to  complex  poles  for  best  cancel  compensation? 
The  search  for  an  answer  to  this  question  initiated  this  study.   Such 
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Effect  of  a  pair  of  Complex  Poles  near  the  Imaginary 
Axis  on  High  Gain  Systems 


%  : 


+  o, 


m ygg— * — 


■zz± 


-wH-t-fc 


mi 


— k 


ooo 


Fig.  2 

Hoot  Loci  from  Complex  Pole  with  adjacent  Complex  Zero 


a  location  or  area  might  make  it  possible  to  locate  complex  zeros  for 
the  design  of  self-adaptive  compensators. 

First  in  this  study,  the  angle  of  emergence  from  a  complex  pole 
and  zero  was  studied  to  determine  how  it  varied  and  how  it  was  related 
to  the  transfer  function. 

Nextj 43  different  case  studies  were  programmed  into  the  Control 
Data  Corporation  1604  Computer  to  obtain  215  root  locus  plots.   These 
plots  and  observations  are  contained  in  Chapter  3. 

The  term  "compensated  system"  will  be  defined  hereafter  to  mean  a 
system  with  a  pair  of  complex  conjugate  zeros  relatively  close  to  com- 
plex conjugate  poles.   The  "uncompensated  system"  will  be  defined  as 
the  system  without  the  complex  conjugate  zeros,  but  will  include  the 
real  poles  associated  with  a  "zero-compensator." 
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CHAPTER  II 
.  EMERGENCE 

1,   Determination  of  the  Angle  of  £;:er,~er.ce 

The  angle  of  emergence  of  the  root  locus  from  a  complex 
pole  or  zero  has  been,   defined  in  the  literature  as  the  angle 
of  the  tangent  to  the  root-locus  curve  at  the  pole  or  zero  as 
shown  in  lrlg.  3»   A  test  point  s-,  near  the  complex  pole  on 


the  root  locus  satisfies  the  ecuation 
N. 
Z.   -  \P, 


where  n  is  an  odd  integer  or  zero,  N  is  the  total  number  of 
zeros,  M  is  the  total      r  of  poles,  and/? ^   is  the  angle 
of  the  complex  pole  infinitesimally  close  to  the  test  point. 
Fig.  3  shows  that  /?-,   is  measured  counter-clockwise  with  the 
negative  abscissa  axis  as  reference,  whereas  the  angle  of 
emergence,  (]/pis  measured  counterclockwise  with  the  positive 
abscissa  axis  as  reference.   It  is  thus  seen  that  the  angle 
of  emergence  is  equal  to  the  angle  from  the  Lest  point  to 
the  complex  pole,  using  their  corresponding  references.  Thus, 

The  last  two  terms  in  the  above  ecuation  can  be  deter- 
mined with  the  complex  pole  or  zero  (from  which  the  angle  of 
emergence  is  to  be  determined)  as  the  origin  of  the  co- 
ordinate system,  since  the  test  point  s^  is  infinitesimally 
close  to  the  complex  pole  (or  zero).   Thus,  the  origin  of  the 
s-plane  is  translated  to  the  desired  complex  pole  or  zero  by 
setting  w  =  s  +  °*— j  &  where  w  is  the  new  complex  variable  in 
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Fig.  3 

Angle  of  Emergence  of  complex  pole,  Of 
Angle  of  Emergence  of  complex  zero,  ((4- 
Angle  from  test  point  a-^_   to  comolex  pole,  L£l 
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the  w- plane.   A  test  point  wn  near  the  new  origin  (pole  or 
zero)  which  is  on  the  root-locus  must  still  satisfy  the 
angle  criterion  (130  degrees  or  an  odd  multiple)  .   As  )w| 
approaches  zero  the  sum  of  the  angles  on  the  root  locus  plus 
180  degrees  would  be  tae  angle  from  the  origin  to  the  test 
point.   When  ]wj is  identically  zero,  there  is  no  angular  con* 
trioution  from  the  pole  or  zero  at  the  origin  and  thus  the 
sum  of  the  remaining  angles  is  the  angle  of  emergence.   For 
example,  a  characteristic  equation  is  partitioned  as  follows 

s(S  +  10Ms  +  2-u>y   "  -    "     , 

To  find  the  angle  of  emergence  at  the  complex  pole  located 
at  —2  +j5  in  the  s-plane,  transfer  the  origin  to  this  pole 
by  setting   w  =  s+2  -j5  or  s  =  w  -  2  +j5    Then, 

(w-2+j5)('w-      """  -  .  """  -:5)(v/-:  ,  _-2-:5) 


(w-2  +  j5)(v/+b-ei5Uv/)Cw+jlO) 
(w+1) (w+1- 

(  -1 )  ( v.-2  -  J,)  { v/+  j  10  )  (  v,-t6  +  j  5  ) 


(2-2) 


Define  the  angle  of  emergence  from  the  comolex  oole, 

<*=  &,w»«0  (2"2a) 

S  =  tan"1  §  +  tan"1  1°   -  (180  +  tan"1,!  +  tan"1  § 
+  90) 

ty?=  0  +  84.5  -  (180  +  111.6  +  32  +  90) 

tf,=  30.9 
Now  to  find  the  angle  of  emergence  from  the  complex  zero 
located  at  -3  +J5  in  the  s-plane,  transfer  the  origin  to 
this  zero  by  setting  w=s+3_J5   or  s=w-3+j5 
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(w-3+.i5  +  3-.i5)(w-3-M5-H+.i5) —  =  _! 

(w-3  +  d5)(w-3+j5+10)(w-3+j5+2-d5)(w-3+d5+2+dt>; 

(w)(w+,jlO) _  =  -i 

(w-3+j5)(w+7+j5)U/-D(w-l+olo; 


(-1)  (w-3+.i5)(w+7  +  .i5)(w-l)(w-l  +  ,jlO) 
w (w+jlO) 


(2-3) 


Define  the  angle  of  emergence  from  the  complex  zero, 


jwj  =  0 


-1  5 


-1  5 


^=  180  +  tan   £  +  tan   y  +  *-«  _x 

-  90 
(£=  180  +  121  +  35.5  +  180  +  95-90 


-10    .   -1  10 
an   v  +  tan   ^y 


(2-3a) 


2.   Complex  Residues  and  the  Angle  of  Emergence 

The  closed  loop  transfer  function  can  be  expanded  by 
partial  fraction  expansion  and  the  residues  evaluated  by  use 
of  Heaviside's  expansion  theorem:  (a)  multiply  function  by 
(s  -  p.  ) ,  and  (b)  let,  s  =  p  ,  giving  the  residue  at  K  .  • 
Using  the  previous  example  again,  the  residue  angle  at 
s  =  -2  +j5  is: 


4   -  (s+.V-.i5)(s+3+,i5) 
\c     s(s+10)(s+2+j5) 


'+J5 


(2-4) 


It  is  seen  that  this  residue  evaluation  is  similar  to  eoua- 
tion  (2-1)  except  for  the  -1  in  the  denominator  which  intro- 
duces 180  degrees.   Thus  the  angle  of  the  residue  at  a  com- 
plex pole  differs  from  the  angle  of  emergence  by  180  degrees 


<i= 


K 


+  180 


(2-5) 


The  characteristic  equation  of  the  preceding  section 
may  be  partitioned  as  follows,  where  the  zeros  are  now  the 
"poles"  and  the  poles  are  now  the  "zeros": 
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The  left  side  of  this  equation  is  seen  to  be  the  reciprocal 
of  the  closed-loop  transfer  function.   Although  this  is  an 
improper  fraction,  the  residues  can  still  be  evaluated  by 
Keavlside'e  expansion  theorem.   Thus,  the  residue  at 
s  =  -3+J5  is: 

Kz  "  (r;  +  3-j5)  s  =  -3+j5         U"6 } 


L 


L 


.    ( -3  +  .15)  (-3^n^]0)  (-3+.15-?-,15  )(-">-M5  +  ?-M5)    (2-7) 
3  (-3-j5+3+:5) 


/. 


K7  =  341.5 


Comparing  equations  (2-7)  and  (2-3)  it  is  seen  that 

they  differ  by  130  degrees.   Thus  the  angle  of  emergence 

from  a  complex  zero  is 

(fe  =  /\  *  l8°  (2-7a) 

where  IlL 


is  determined  as  illustrated  above. 


3-  Effect  of  Comjlex  Zeros  near  Complex  Poles 

What  effect  does  the  relative  position  of  complex  zeros 
adjacent  to  complex  poles  have  on  the  angle  of  emergence 
from  the  complex  poles?  To  answer  this  question, first  con- 
sider the  change  of  the  angle  of  emergence  with  the  zeros  at 
constant  racius  from  the  poles  while  varying  the  angular 
oosition  of  the  zeros  with  resoect  to  the  complex  poles. 

It  is  assumed  that  the  displacement  between  complex 
poles  and  the  distance  between  complex  pole  and  its 
neighboring  zero  is  of  a  ratio  of  approximately  10:1.  As 
seen  in  Fig.  4,  as  the  angle  9,  the  angular  position  of  the 
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Fig.  4 

2frect.  of  Zero  Angular  Position  on  Angle  of  Emergend 


16 


complex  zero,  increases,  the  angle  of  emergence  increases 
in  magnitude.   The  angular  contribution  of  X  due  to  the 
lower  half  plane  zero  is  negligible  compared  with  G  in  the 
upper  half  plane.   Therefore  a  linear  slope  of  one  exists 
between  angular  position  of  zero  and  angle  of  emergence,  or 

3e  (2-3) 

Next  consider  Fig.  5  which  snows  the  effect  of  radial 
position  of  complex  zeros  (with  constant  angular  oositlon) 
on  the  angle  of  emergence.   As  the  zeros  vary  in  radial  posi- 
tion, there  is  no  angular  change  cue  to  upper  half  plane  zero. 
The  change  in  lower  naif  plane  zero  angle  is  0.   In  addition 
to  the  assumption  stated  ab  ve,  the  line  P^Zp  is  assumed 
perpendicular  to  line  PpA  for  small  changes  in  radial  posi- 
tion r.   Then 

0   =  tan"1/  r  cos  Q    \  where   d  is  the  displacement 

[&  -   r  sin  Qj     distance  between  complex  poles, 


pl?2 


=ife 


2     2 

1  +  r   cos  &    p 

(d  -  r  sir  G~) 


1£-    ..    f(r   cos    Q)(-l)Ql-r    rir.   GN         -    :  -,    G)    +  (d-r   sinGl^cosQ] 


y   - 


-  2        2r 

1    +        r   cos   G 


[d    -   r   sin  W) 

-*.u        fr    sinG   cos    G    _      +  cos   G  "? 

B"=    L,-    -    V    ^.r.   «>2  d    -   r   s^    6    / 

2        2 

1      +        r  cos   6-  p 

(d   -r   sin   G) 
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Fis.  5 

Effect  of  Zero  Radial  Displacement  on  Angle  of  Emergence 
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~a<f>  r     .  ;  r    G    o  o  r,    G      +      o  o  s    6  (  d    -   r    si  r.    G ) 

?/    "    (d    -   r   sin   0)^      +      r2cos20 

^    d  nor.  G 

"^  ~~  ~" ?  I     2 

d   -  2dr  sin  9   +  r 

2  2 

Sine  it  was  assumed  that  d/r  =  10,  d  >>  -2drsin9  +  r 

>£  _   :os  Q  ^2~9) 

-^  "     d 

The  total  derivative  of  the  angle  of  emergence  is: 

For  smal 1  d i  s pi ac ement s , 

(^  =  0   +  v<;',,    ft    +  ^  (2-10) 

d 

compensated 
/Si  it .   If  r, 

adjacent  complex  poles  and  zeros,  is  allowed  to  oecoxe  zero, 

or  "perfect  cancellation"  with  poles,  the  anguler  position 

of  the  zeros  is  undefinea.   The  constant  ^is  then  evaluated 

as  the  sua;  of  all  remaining  zeros  in  the  system  minus  the 

sum  of  all  remaining  poles,  measured  from  the  upper  half 

plane  "cancelled"  complex  pole.   Thus  in  the  example  cited 

in  section  1,  A   is  evaluated  as  210.9.   For  9  =  180,  the 

angle  of  emergence  04,  =  30.9  degrees •   The  second  term  on  the 

right  side  of  equation  (2-10)  is  negligible  compared  with 

the  other  terms  under  the  assumptions  that  have  been  made. 

(See  Fie.  7)   It  is  seen  that  the  angled  differs  from  the 

residue  angle  of  the  uncompensated  system  (no  complex  zeros 

present)  by  90  degrees. 

K   =  A-   90  (2-11) 

p 
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The  angle  of  emergence  of  the  uncompensated  system, 
(j/  ,    is  readily  determined: 

^  =  K   +  180  (2-5a) 

0/  =  /2      -   90  +  180  =/&   +  90  (2-12) 

Equation  (2-7)  can  now  be  written  as: 


(£  =  G 


cosO 


V       + 


l/p  -  90 


L/c    S  Q   +  VP  -  90  (2-13) 

Thus  the  angle  of  emergence  changes  directly  with  the 

zero  position  angle  0.  (See  Fig.  8) 

The  residue  angles  of  the  complex  pole  and  zero  are 

defined  in  equations  (2-4  )  and  (2-7)  respectively  and  can 

be  written  in  the  following  form:: 

4 


A 


L £—         ~/C 


at  pole  p. 


it  zero  z, 


(2-14) 
(2-15) 


per- 


The  relationship  between  the  pole  residue  angle,  K 

and  the  zero  residue  angle  K  can  now  be  written  as: 

z 

.A  =  A=   +Zfi  +Z£^  (2_16) 

Letx=    /z0      raid'  ^  =    /n  The    equation    (2-16)    becomes: 

A* =  -Aoc  +  •*  +  7  (2-17) 

The  angles  •*  and  ^  are  seen  in  Pig. 6.   It  is  observed  from 
this  figure  that  <*  =  G  -  180   and  ^  =  G 

Substituting  these  relationships  into  equation  (2-17) 


yields  the  following  result: 


EC   =   -A 
z     // pj 


+  2G  -  180 


(2-18) 
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Fi£.  6 


Relationship  between  Residue  Angle  of  Complex  ?ol< 
ana  Complex  Zero. 


/4  '-    4k 


cX      + 


f 
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But  from  equations  (2-5)  and  (2-7a), 


L 


K   =  (/c  -   180 
pc 


A7  =  <&  -   180 

Therefore, 

(^  _   180   =  -   U/c  +  180   +   2G  -   180  or, 

^  =  _  </£  +   2G   +  180  (2-19) 

Substituting  equation  (2-13)  into  equation  (2-19)  yields 
the  relation  between  the  angle  of  emergence  from  the  com- 
plex zero  and  the  angle  of  emergence  from  the  complex  pole 
when  the  complex  zeros  are  not  present. 

Q£=   -9  -  <^+  90  +  2G  +  180 

(f£=     Q     -  Of-   90  (2-20) 

Using  the  relation  from  equation  (2-5a),  then 

££=  9  - /\     +   90  (2-21) 

Subtracting'  equation  (2—13)  from  equation  (2-20)  gives 
the  relation: 

(£-&=  -  ?(/?=   -2(/frp  +   180) 

t^-</c   =   -2/^  (2-22) 

Thus  the  difference  between  the  angle  of  emergence 
from  a  complex  pole  and  the  angle  of  emergence  of  its 
neighboring  complex  zero  is  a  constant  for  a  given  system, 
and  thus  not  a  function  of  zero  position  0. 

Figures  7  and  8  are  computer-solution  graphs  of  the 
angle  of  emergence  versus  complex  zero  radial  and  angular 
positions  respectively* 
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CHAPTER  III 
ROOT  LOCI  LEAVING  COMPLEX  POLES 

1.    Data  Designation 

Having  analyzed  the  angles  of  emergence  from  complex  poles  and 
zeros,  the  next  step  in  this  study  is  to  determine  the  effect  of 
various  system  transfer  functions  on  the  root  loci  emerging  from  the 
complex  poles,  and  to  determine  when  this  loci  closes  on  its  adjacent 
complex  zero,  i.e.,  the  root  loci  starts  at  the  complex  pole  and  ends 
on  the  complex  zero.   As  shown  in  Fig. 2  of  Chapter  I,  this  loci  may 
go  to  the  zero  at  infinity  instead  of  the  complex  zero. 

Over  fifty  different  cases  were  tried  and  the  data  consisted  of 
about  300  root  locus  plots  using  the  computer  program  in  Appendix  I. 
Only  43  case  studies  (Fig.  9)  are  presented  with  215  individual  root 
locus  plots  (Fig.  10)  for  this  analysis. 

Each  case  study  or  individual  root  locus  plot  is  designated  or 
classified  using  letters  and  numbers,  such  as  "A-la",  where  the  first 
(capital)  letter  designates  generally  the  number  of  real  poles  and 
zeros,  and  position  of  complex  conjugate  poles.   The  number  indicates 
the  variable  position  of  the  given  number  of  real  poles  and  zeros. 
The  last  small  letter  specifies  the  variable  position  of  the  complex 
zero  with  respect  to  the  complex  pole. 

Double,  triple,  etc.,  real  poles  are  specified  by  placing  an  "X" 
over  an  "X"  or  a  group  of  "X's"  near  appropriate  point  on  the  plot. 
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The  double  real  poles  are  arbitrarily  used  to  represent  the  poles 
associated  with  a  "cancel-compensator"  and  is  usually  located,  for 
this  study,  at  -8  units  except  in  a  few  cases  they  are  placed  at  -4 
un  i  r  a . 

Case  A  is  a  fifth  order  system.   Case  B  is  a  fifth  order  system 
with  one  real  zero.   Case  C  is  a  sixth  order  system.   Case  D  is  a 
sixth  order  system  with  the  complex  conjugate  poles  moved  out  to  -3 
+  j4.   Case  E  is  a  sixth  order  system  with  one  real  zero.   Case  F  is 
a  seventh  order  system.   Case  G  is  a  seventh  order  system  with  one 
real  zero.   Case  H  is  also  a  seventh  order  system  with  one  real  zero. 
Case  I  is  an  eighth  order  system.   Figure  9  presents  these  43  case 
studies,  and  shows  only  the  pole-zero  configurations  without  the  com- 
plex conjugate  zeros.   Figure  10  shows  the  root  locus  plots  for  these 
case  studies  when  the  complex  zeros  are  placed  in  various  positions. 

Table  1  is  a  summary  of  the  graphical  data  from  Fig.  10.   In 
addition  it  contains  the  computed  values  of  angles  of  emergence  from 
pole  and  zero  (for  9  s  0°)  and  the  uncompensated  residue  angle  using 
the  formulas  developed  in  Chapter  II.  A  comparison  with  these  calcu- 
lated values  from  the  transfer  function  compares  very  good  with  com- 
puter root  locus  plots. 
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2.    Observations  from  Data 

a.  There  is  always  a  zero  position  angle  9  in  the  first 
quadrant  with  respect  to  the  complex  pole  which  will  cause  the  root 
locus  from  the  complex  pole  to  end  on  the  near  complex  zero.   This  is 
observed  in  all  cases. 

b.  If  the  root  locus  closes  on  the  complex  zero,  then  the 
locus  from  the  pole  moves  in  direction  of  least  number  of  degrees 
(angle  of  emergence  from  pole  compared  with  angle  of  emergence  from 
zero) . 

c.  The  locus  closes  on  complex  zero  if  there  is  a  root  locus 
segment  on  negative  real  axis  to  the  right  (closest  to  the  origin)  of 
the  real  part  of  the  complex  conjugate  pole.  (C-l,  C-2,  D-l,  D-2,  D-3) 

d.  Just  as  the  radial  distance  of  zero  had  little  effect  on 
angle  of  emergence,  it  also  has  little  effect  on  root  locus  plots  ex- 
cept to  make  plots  proportionately  larger  or  smaller.   (Compare  Case 
A-l,  A-5,  and  A-6) . 

e.  The  locus  from  complex  pole  closes  on  complex  zero  if 
asymptote  vector  is  90  degrees  and  the  asymptote  centroid  is  far 
enough  out  on  negative  real  axis,  (i.e.,  Number  of  pole  minus  number 
of  zeros  is  equal  to  two).   See  cases  B-l,  B-2,  B-3,  and  B-4. 
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TABLE    1 
SUMMARY   OF  DATA* 


Dcsig 


;C'  Transfer  Function^    Other  Data 
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*"Desig"  is  the  designation  of  each  root  locus  plot,  0  (theta)  is  the 
angular  position  of  the  complex  zero  measured  with  respect  to  its 
adjacent  complex  pole,  "C"  indicates  that  the  root  locus  from  the 
complex  pole  closes  or  ends  on  its  adjacent  complex  zero,  "NC"  in- 
dicates that  the  root  locus  does  not  close  on  the  complex  zero.  Values 
of  /k  were  found  by  use  of  a  spirule,  and  the  values  of  ^  and  (^ 
were  cSlculated  by  use  of  the  formulas  developed  in  Chapter  II  with  9 
equal  to  000  degrees. 
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Table   1    (Cont'd) 
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Transfer  Function:  Other  Data 
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Table  1  (Cont'd) 
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Transfer  Function;  Other  Data 
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Table  1  (Cont'd) 
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Transfer  Function:  Other  Data 
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F-2e 
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Table  1  (Cont'd) 
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Transfer  Function:  Other  Data 
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CHAPTER  IV 
CONCLUSIONS 

Two  analytical  methods  for  obtaining  the  angle  of  emergence  from 
a  complex  pole  or  zero  have  been  presented.   In  the  first  of  these 
methods,  the  origin  was  transferred  to  the  complex  pole  or  zero  (w- 
plane)  and  then  the  angle  of  emergence  was  defined  by  equations  (2- 
2a)  and  (2-3a). 

The  second  method  was  to  determine  the  residue  angle  by  means  of 
Heaviside's  expansion  theorem  and  then  simply  add  180  degrees  to  the 
residue  angle  to  obtain  the  angle  of  emergence. 

The  following  equations  are  exact,  i.e.,  no  approximations  were 
used. 

<£    --/j^.     +    IS°  (2-5) 

<l£   =/Vg   ■+  tSO  (2-7a) 

The  easiest  method  to  obtain  the  angle  of  emergence  is  to  find 
the  residue  angle  with  a  spirule,  then  add  180  degrees. 
Other  relations  developed  in  Chapter  2  are: 

(^   r  O    -hVf-JO  (2-13) 

(£   =  -  (/c  +   ^e  +1VO  (2-19) 

(£-<^    =    -Jz/fr  (2-22) 

and   since  C^p    =  /Kp   ~  /  &0 

Vc    =     S    +Af      *    <?0  (3-D 
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It  must  be  remembered  that  these  equations  are  only  approximate 
with  the  linearizing  assumptions  and  approximations  which  have  been 
made.  Figures  7  and  8  are  computer  solutions  of  the  angle  of  emergence 
vs.  radial  position  of  complex  zero  and  the  angle  of  emergence  vs.  zero 
poi'itionangle  respectively,  and  thus  verifies  the  linear  approximations 
for  the  assumed  separation  distance  between  complex  conjugate  poles. 

Thus  knowing  the  angle  of  emergence  from  a  complex  pole  (^  ,  the 
position  angle  0  of  the  neighboring  complex  zero,  the  angle  of  emer- 
gence from  the  complex  zero  is  readily  obtained  by  equation  (2-19). 
Similarly  if  the  angle  of  emergence  from  a  complex  pole  of  an  un- 
compensated system  is  known,  the  effect  of  adding  complex  zeros  on  the 
angle  of  emergence  is  obtained  by  using  equation  (2-13).   Equation  (2- 
22)  reveals  that  the  difference  between  the  angles  of  emergence  from 
complex  pole  and  zero  is  proportional  to  the  "uncompensated"  residue 
angle  and  is  independent  of  zero  position  angle  9. 

The  angle  of  emergence  from  complex  pole  and  zero  help  to  pre- 
dict the  manner  in  which  the  loci  leave  the  pole  or  zero.  A  vector 
is  drawn  from  the  pole  in  direction  of  the  angle  of  emergence  and  a 
second  vector  is  drawn  from  the  zero  in  the  direction  of  its  angle 
of  emergence.  Then  the  tails  of  these  vectors  are  extended  until 
they  intersect  at  a  point  and  thus  defines  an  angle  less  than  or  equal 
to  180  degrees.  It  was  observed  that  whenever  the  root  locus  from  a 
pole  closes  on  its  adjacent  zero9  the  loci  moves  in  such  a  direction 
so  as  to  be  contained  within  the  angle  defined  by  the  angle  of  emer- 
gence vectors. 
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It  appears  that  there  is  always  a  zero  position  angle  between  0 
and  90  degrees  (the  first  quadrant  with  respect  to  the  complex  pole) 
such  that  the  root  locus  from  the  complex  pole  closes  on  the  complex 
zero.  Thus  a  self-adaptive  zero  compensator  could  be  located  within 
this  90  degree  sector  and  effectively  "compensate"  the  poles.  The 
residues  from  the  roots  of  these  closed  root  loci  would  be  negligible 
compared  with  the  dominant  roots  of  the  control  system. 
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APPENDIX  I 

Two  digital  computer  programs  were  used.   The  first  was  "Program 
ANQ1ERG"  to  determine  the  angle  of  emergence  for  various  pole  and 
zero  locations.   Figures  7  and  8  were  outputs  of  this  program  which 
verified  analytical  developments  of  Chapter  II. 

The  second  digital  computer  program  was  "Program  Eqn"  which  is 
a  program  designed  to  plot  root  locus  for  any  given  characteristic 
equation  (up  to  order  30).   This  program  was  originally  programmed  by 
LT  Joseph  D.  Fennick,  Jr.,  USN.   A  few  changes  in  the  original  sub- 
routine were  made  to  facilitate  data  input,  reduce  printout,  and  speed 
up  graph  plotting  time.   Figures  11,  12,  and  13  are  sample  outputs  of 
this  program.   Table  2  contains  all  of  the  input  data  (characteristic 
equations)  for  each  of  the  215  individual  cases  studied. 
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.. JOB0414F,REES,  E.G. 

PROGRAM  ANGMERG 
C      THIS  PROGRAM  CALCULATES  THE  ANGLE  OF  EMERGENCE  FOR  VARIOUS 
C      POLE  AND  ZERO  LOCATIONS.   MULTIPLE  RUNS  CAN  BE  MADE  BY 
C       INDICATING  NUMBFR  OF  RUNS  AND  ADDING  ADDITIONAL  DATA  CARDS 
C      RFAL  AND  IMAGINERY  PARTS  OF  POLFS  AND  ZEROS  ARE  REAH  IN 
C      WITH  STIPULATION  THAT  REP(l)  AND  YMP(l)  BE  COMPLEX  POLE  ABOUT 
C      WHICH  THE  ANGLE  OF  EMERGENCE  IS  CALCULATED.   COMPENSATION 
C      ZEROS  ARE  NOT  READ  IN*  BUT  ARE  GENERATED  IN  THE  PROGRAM. 
C      LET  REP(2)  AND  YMP ( 2 )  BE  CONJUGATE  OF  RFP(l)  AND  YMP(l) 
C      PART  OF  THE  PRINT  OUT  AND  GRAPHS  CAN  BE  SUPPRESSED  BY 
C      CHANGING  VALUES  OF  W  AND  V  RESPECTIVELY.   FOR  A  CURVE  OF 
C      RADIUS  VS  ANGLE  OF  EMERGENCE  AT  A  PARTICULAR  ANGLE  OF 
C      COMPENSATION  ZERO  FROM  COMPLEX  POLE*  ANGLE  AAG  IS  READ  IN 
C       IN  RADIANS.   TWO  GRAPHS  OF  ANGLE  OF  EMERGENCE  VS  ANGLE 
C      POSITION  OF  COMPENSATING  ZEROS  AT  A  GIVEN  RADIUS  IS  AVAILABLE. 
C      RR1  AND  RR2  ARE  THE  RADII  FOR  THESE  CURVES.  ADDITIONAL 
C      SYMBOLS  ARE —  CE  TROID,  THE  INTERSECTION  OF  ASYMPTOTES*, 
C      THETA  —  ANGLE  OF  ASYMPTOTES  FROM  HORIZONTAL  WHEN  COM  PENSATION 
C      ZEROS  ARE  ADDED* ♦  ZETA  -  ANGLE  OF  ASYMPTOTES  FROM  HORIZONTAL 
C      WITHOUT  COMPENSATING  ZEROS**  ANGZ  -  SUMMATION  OF  ANGLES  TO 
C       ALL  ZEROS  FROM  COMPLEX  POLE,,  ANGP  -  SUMMATION  OF  ANGLES  TO 
C      ALL  POLES  FROM  COMPLEX  POLE. 

C      DATA  READ  IN  AS  FOLLOWS 

C  NOS     12     NUMBER  OF  RUNS 

C  M,N     212      NO.  OF  POLFS  M,  NO.  OF  NON.  COMP.  ZEROS  N 

C  REP(I)    8F10.3     REAL  PART  OF  POLES 

C  YMP(I)    8F10.3       IMAG.  PART  OF  POLES 

C         REZ(I)     8F10.3     REAL  PART  OF  ZEROS 

C  YMZ(I)        8F10.3     IMAG.  PART  OF  ZEROS 

C  RR1,RR2,  AAG,  W,  V     5F10.3   RRl  AND  RR2  DESIRED  RADIUS  FOR 

C  •        GRAPH,  AAG  IS  ANGLE  IN  RADIANS  LOCATING  POSITION  OF 

C  COMPENSATING  ZERO  FOR  RADIUS  VS  ANGLE  OF  EMERGENCE  GRAPH 

C  W  SUPPRESSES  PRINTOUT  AS  FOLLWS  (W  =  0.0  PRINTS  ONLY  EIGHT 

C  VALUES,  W  =  1.0  PRINTS  OUT  SIXTEEN  VALUES  FOR  EACH  RADIUS) 

C  V  SUPPRESSES  GRAPHS  AS  FOLLOWS  (  V  =  0.0  ALL  GRAPHS 

C         SUPPRESSED,   V  =  1.0  A  TOTAL  OF  THREE  GRAPHS  PER  RUN 

DIMENSION  RFP( 8) ,YMP( 8  )  ,REZ(8 ) ,YMZ (8 ) ,R(900) ,B(900) ,REZ1 (900 )    » 
1ZANG(900) ,ANG(900) • I  TITLE! 12) ,YMZ1 (900) ,ZA(8) ,PA(8)  ,GNA(900)  , 
2  D(900),  AZ(QnO),  CEN(900),  ACF(900) 

READ  63*  NOS 
63  FHPMAT(  12) 

DO  402  K=  l,NOS 

READ  2,M,N 
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FORMAT 

(  212  ) 

READ  4, 

(REP( I ) ♦ I  = 

1,M) 

READ  4, 

( YMP( I ) , 1= 

1  ,M) 

IF(N-1  ) 

4,9,9 

READ  4* 

(REZ(  I  )  ,  1  =  1 

,N) 

RFAD  4* 

( YMZ(  I  )  ,  1  = 

1,N) 

FORMAT 

(8F10.3) 

READ  44 

,  RRl,  RR2, 

AAG,W,V 

FORMAT 

(5F10.3)   * 

PRINT  DATA  READ  IN 

PRINT  63»K 
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62  FORMAT  (/4X,7HRUN  NO. ,12) 

PPINT  6»M,N 
6  FORMAT ( //9X,2HM=,I3,20X,2HN=,I3) 

PRINT  8 
8  FORMAT  (/5X.7H  REP(I),9X,7H  YMP(I),9X,7H  REZ(I),9X,7H  YMZ (  I  )  ) 

IF(N-l)  48,49,49 

48  PRINT  99,  ( REP (  I  )  , YMP (  I  )  ,  1=  ],M) 
GO    TO  99 

49  PPINT  10,  (RFP(  I  )  ,YMP(  I  )  ,REZ (  I  )  ,YMZ(  I  )  ,  1  =  1  ,N) 
LL  =  N  +  1 
PRINT  99,  (REP(I),  YMP ( I ) ,  I=(_L»M) 

99  FORMAT  (2F15.3) 
10  FORMAT  (4F15.3) 

PRINT  45,  RR1 ,RR2,AAG,W,V 
45  FORMAT  (/4X,5HRR1  =,F10.3»  5X,5HRR2  =,F10.3,  5X,5HAAG  =,F10.3, 
15X,  3HW  =,F10.3,5X,3HV  =,F10.3) 

XX  =  M-N 

X  =  XX  -  ?. 

ASYMPTOTE  VECTOR  WITH  COMPENSATING  ZEROS,  THETA 

IF  (X)  202,202,203 

202  PRINT  205 
20^  FORMAT( /2X, 1 SHTHFTA=  INFINITY) 

GO  TO  2  06 

203  THFTA  =  180. /X 
PRINT  204,  THETA 

204  FORMAT  ( / 2X ♦ 6HTHF T A=  ,  F10.3) 
ASYPTOTE  VECTOR  WITHOUT  COMPENSATING  ZEROS,  ZETA 

206  I-F  (M-N)  210,  210,  211 

210  PRINT  212 
212  FORMAT  ( /2X , 14HZET A=  INFINITY) 

GO  TO  240 

211  ZETA  =  180. /XX 
PRINT  214,  ZrTA 

214  FORMAT  ( / 2X , 5HZ ETA= , F 1 0  .  3 ) 
C      ASYMPTOTF  CENTROIO,  CEN  =  (SUM  -POLES  -  SUM  ZFROS)/  M-N 
240  SUM=0.0 

DO  201  1=1,8 
201  SUM  =  SUM  +  REP(  I  )-REZ (  I  ) 

SMM=-SUM 
C      ANGLES  DUE  TO  NON-COMP.  ZEROS,  ANGZ 

ANGZ=0.0 

IF  (N)  250,250,  251 
25  1  DO  12  I  =  1,N 

IF(YMZ(  I  )-YMP(  ]  )  )71  ,72,73 

71  •IF(REZ( I )-REP(l ) )74,76,77 

74  ZA(I)=ATANF(  (YMP( 1  )-YMZ(  I  )  )  / ( RFP (  1  ) -REZ (  I  )  )  ) +3. 14159 
GO  TO  12 

72  IF(REZ(  I  )-REP(l  )  )75,80,80 

75  ZA( I )=3. 14159 
GO  TO  12 

76  ZA(  I  )  =  3./2.*3. 14159 
GO  TO  12 

7  7  ZA(  I  )=2.*-3.  1415  9-ATANF(  (YMP(  1)-YMZ(  I  )  )/  (REZ(  I  )-REP(  1  )  )  ) 

GO  TO  12 
80  ZA(  I  )=0.0 

sa  te  ia 
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73  IF(REZ( I )-REP(l ) )81,82,83 

81  ZA(I)=ATANF(  ( YMZ(  I  ) -YMP ( 1 )  )  / ( RFZ (  I  ) -REP ( 1 )  )  ) 
GO  TO  12 

82  ZA (  I  )  =  3. 14159/2. 
GO  TO  12 

83  ZA(I)  =  ATANF(  ( YMZ (  I  )-YMP( 1 )  )  / ( REZ (  I )  -  REP(l))) 
12  ANGZ=ANGZ+ZA( I ) 

ANGLES  DUE  TO  POLES,  ANGP 
250  ANPG  =  0.0 

IF  (M-2)  254,  254,  253 
253  00  14  I  =  3,M 

IF( YMP(  I  )-YMP( 1)  ) 9 1,9 2,  9 3 

91  IF(REP( I )-RFP(l ) )94, 96,97' 
9  4  PA( I )=ATANF( (YMP( 1)-YMP( I )  ) / ( RFP ( 1 ) -REP ( I )  )  ) +3 .  141 59 

GO  TO  14 

92  IF(REP( I )-REP(l ) )95, 100,100 

95  PA(  I  )  =  3. 14159 
GO  TO  14 

96  PA(  I  )=3./2.*3. 14159 
GO  TO  14 

97  PA(I)=2.*3. 14159-ATANF(  ( YMP ( 1 ) -YMP (  I )  ) / ( REP (  I  ) -REP (  1  )  )  ) 
GO  TO  14 

100  PA (  I  )=0.0 
GO  TO  14 

93  IF(REP( I )-RFP(l )) 101 ,102,103 


101  PA(I)  =  3.14159  -  ATANF( ( YMP( I )-YMP( 1 ) )/ (REP( 1 )  -  REP(I))) 
GO  TO  14 

102  PA(  I  )  =  3. 14159/2. 
GO  TO  14 

103  PA(I)  =  ATANF(  ( YMP (  I  )  -YMP ( 1 )  ) / ( REP (  I  ) -REP ( 1 )  )  ) 

14  ANPG  =  ANPG  +  PA( I ) 
ANGP  =  ANPG  +  3.14159/2. 
GO  TO  2  52 

254  ANGP  =  3.14159/2. 
C      ANGLE  OF  EMERG.  IF  COMPLEX  ZEROS  NOT  PRESENT,  ANGLE 

25  2  ANGLE=180./3.14159*(-3.14159+ANGP-ANGZ) 
C      GENERATE  COMPENSATION  ZEROS 
C      ALT.  STATE.  15-19  FOR  VARIABLE  RADIUS 
DO  301  J=l,20 
CA  =  J 
R( J)=CA/10. 

15  DO  26  1=1,630 

16  AI=I-1 

17  B(  I  )=AI/100. 

18  RFZ1(I)=  REP(l)  -  COSF(B( I ) ) *R ( J) 

19  YMZKI)  =  YMP(l)  +  SINF(B(  I  )  )*R(  J) 

C      ANGLES  DUE  TO  COMPENSATION  ZEROS,  ZANG(I) 

53  IF(B( I )-3. 14159)     54,54,55 

54  D(I)  =  180.    -    B(  I  )*180./3. 14159 
GO    TO    5  6 

55  'D( I )=540.-B( I )*180./3.14159  " 
5  6     IFfREZK  I  )-REP(  1)  )  5  7,5  8,59 

5  7  "ZANG(  I  )=180./3.14159*(3.14159-ATANF(  (-YMZK  I  )-YMP(l)  )/ 
1(REP(  D-REZK  I  )  )  )  )+D(  I  ) 
GO    TO    20 
&8   -2AN6(  I  )=96.*ai  I  \ 
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GO  TO  20 
5  9  7ANG(  I  ) =180./3.14159*( ATANF(  (-YMZ1 (  I )-YMP( 1 )  )/(REZl  (  I  )-REP(  1)  )  )  ) 
•1  +  0  (  T  ) 

r  of  EMERGENCE  WITH  COMPENSATION,  ANG(I) 

20  GNA (  I  )=ANGLF-ZANG(  I  ) 
IF  (X)  229,229,230 

230  rcN(I)  =  (SMM  +2.  *  RFZKI)  )/X 

229  AZ(I)  =  180.  /  3.14159  *  (6.28  -  BID) 

21  IF(GNA( I ) )  22,23,23 

22  ANG(  I  )  =  -GNA(  I  ) 

IF  (ANG(I))  600,333,333 
600  ANG( I )  =  360.  +  ANG( I ) 
GO  TO  3  33 
?^  ANG( I )  =  360.  -  GNA( I ) 
333  IF  (B( I )  -AAG)  26,46,26 
46  ACF{ J)=  ANG( I ) 

26  CONTINUE 
IF(R(J)-RR1)  321,27,320 

320  IF(R( J)-RR2)321,27,321 

27  LA=4HANG 
IF(V)  31,29,31 

29  GO  TO  321 
31  DO  1  1=1,12 
1  ITITLEt I)=8H 

ITJTLF( 1 ) =8H  REFS 

ITITLE(2)=8HE.G. 

ITITLF( 3)=8HANGLE 

ITITLF(4)=8H  OF 

ITITLF(5)=8H  EMFRG 

ITITLF(6)=  8H  VS  ZERO 

ITITLF(7)  =  8H  POSIT 

IF(R(J)-RR1)  330,331,330 

331  ITITLE( 8)  =  8H  R=.5 
GO  TO  332 

330  ITITLF(8)  =  8H  R=l  .0 

332  ITITLEt 12)=  8H  THESIS 

CALL  DRAW  (  630 , A? , ANG , 0 , 0 , L A , I T I TL E , 0 , 0 , 2 , 1 ♦ 
1  0,0,  7,9,  0,  LAST) 

321  PRINT  405,  R(J) 

405  FORMAT  ( / 2X , 8HR AD  1 1 J^  =,  F10.3) 
22  0  PRINT  2  8 

28  FORMAT  (8X,14H  ANGLE  TO  ZERO,5X,  19H  ANGLE  OF  EMERGENCE,  5X, 
1  9H  CENTROID) 
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IF(W) 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRTNT 
PRINT 
PRINT 
PRINT 
PRINT 
PrtiNf 


60, 
30, 
30, 
30, 
30, 
30, 
30, 
30, 
30, 
30* 
30, 
30, 

3D, 


61,60 
AZ( 1 )  ,ANG(  1 ) 
AZ(39 ) 
A7 ( 78)  , 
AZ( 117)  , 
AZ( 157)  , 
AZ( 196)  , 
AZ(235)  , 
AZ(276)  , 
AZ(314)  , 
AZ(353)  , 
AZ  (392) 


CEN( 1 ) 

ANG(39),  CEN(39) 
ANG  (78)  ♦  CEN(78) 
ANG( 117  )  ♦  C^Nt 117  ) 

ANG  (157),  CEN  (157) 
ANGU96)  ♦  CFN( 196) 
AN'r-  (2^»5  )  ♦"  r^N(  2  35  ) 
ANG(276),  CFN(276) 
ANG(314),  CEN(314) 
ANG(353),  C^N(353) 

ANG(392),  CEN(392) 


«^*^H  ANG(431)|  CPNC43U 
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PRINT 

30, 

AZ(471),  ANG(471),  CFN(471) 

PRINT 

30, 

AZ(510),  ANG(510),  CFNI510) 

PRINT 

30, 

AZ(549),  ANG  (549),  fFN(549) 

PRINT 

30, 

AZ(588),  ANG(588),  CFN(588) 

GO  TO 

30 

61 

PPINT 

30, 

AZ( 1) ,ANG( 1 ) ,  CEN( 1 ) 

PRINT 

30, 

AZ(78)  ,  ANG  (78)  ♦  CEN(78) 

PRINT 

30, 

AZ(157),   ANG  (157),  CEN  (157)" 

PRINT 

30, 

AZ(235),  ANG(235),  CFN(235) 

PRINT 

30, 

AZ(314),  ANG(314),  CFN(314) 

PRINT 

30, 

AZ  (392) ,  ANG(392) ,  CEN(392) 

PRINT 

30, 

AZ(471),  ANG(471),  CEN(471) 

PRINT 

30, 

AZ(549),  ANG  (549),  CEN(549) 

30 

FORMAT  (3F20.3) 

301 

CONTINUE 

PRINT 

401 

>  ANGZ,  ANGP 

401 

FORMAT(/6HANGZ  = ♦ Fl 0 . 3 » 5X , 6HAMGP  =»F10.3) 

PRINT 

32, 

ANGLE 

32 

FORMAT  (/26HANGLE  OF  EMERG.  W/O  COMP  =,F10.3) 

C 

IF  V  = 

=  0.0,  NO  GRAPH.   IF  V  =  1.0   A  GRAPH  OF 

IF(V) 

66,65,66 

65 

GO  TO 

402 

66 

LA=4HANG 
LA  =4HANG 

DO  47 

I  = 

1,12 

47 

ITITLE  (I 

=  8H 

ITITLF(  1)  = 

=8H  REFS 

ITITLF( 2)= 

=8HE.G. 

ITITLE(3)= 

=  8HRADIUS 

ITITLE(4)= 

=  8H  VS 

ITITLE(5) 

=  8H  ANGLE 

ITITLE(6) 

=  8H0F  EMERG 

CALL  DRAW 

(20,R,  ACE,0,0,LA, ITITLE, 0,0, 2,1* 

10,0,7 

►  9,0 

>LAST) 

402 

CONTINUE 

67 

END 

1 
4 

END 

1.0 

1.0 

12.0 

-6.0 

6.0 

1.2 

2.4 

-3.77                1.0 

R  VS  ANG  OF  EMERG 
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.J080414F,  RPFS»F.  G. 
PROGRAM 

READ  240,  NUMB 
pun     FOP MAT  (12) 

DO  241   I  =  It  NUMB 
241  CALL  SPLANE 
FND 

ROUTINE  SPLANE 
DIMENSION  R(129),X(1?9),IT(10)  ,ROOTR( 128)  » ROOT  I  (  128  )  • IT  I TLE( 12)  • 
1A(129),R(129) »ROOTJ(  128)  ,ROOTM(  1  28 )  ♦  AP ( 1 29 )  ♦ AZ ( 1 29 ) 
COMMON  P,VAR,NO,ROOTR,ROOTl ,AP,X 

206  MOD=l 
LAR=4H 

68  FORMAT(8E15.5) 

70  FORMAT    ( / , 14H I  MAG  I  NARY  PART,//) 

69  FORMAT     (/,9HREAL  PART,//) 

200  FORMAT(6A8) 

203  FORMAT(  13  ) 

204  FOPMAT(Fin. 5) 

PFAD  200, ( ITITLE( I )• 1=1,6) 
PPINT  324 
324  FORMATf 1H1 ♦ 11HGRAPH  TITLE*/) 
PRINT  200, ( ITITLFf I ) ,1=1,6) 
28  FORMAT   (/,36HORDER  OF  THE  CHARACTERISTIC  EQUATION,/) 
PRINT  28 
READ  203, NO 
PRINT  203, NO 
N=N0+1 

205  FORMAT  (8F  in.',  ) 

207  FORMAT  (8E12.5) 

22  FORMAT   (  /  ,  53HCONST ANT  PART  OF  THE  COEFFICIENTS  IN  DESCENDING  ORDE 
1  R  ,  /  ) 

PRINT  22 

RFAD  205, ( A(K) ,K=1  ,N) 

PRINT  207, ( A(K) ,K=1 ,N) 

23  FORMAT   ( /,53HVARIABLE  PART  OF  THE  COEFFICIENTS  IN  DESCENDING  ORDE 
1R,/) 

PRINT  23 

READ  205, (B(K) ,K=1 ,N) 

PRINT  207,  (B(K)  ,K=1 »N) 

VAR  =  .001 

ND  =  8 

X^CALF  =  2.0 

yc,CALF  =  XSCALE 

201  FORMAT  (/////, 21H0THE  SYSTEM  POLES  ARE,/) 
PRINT  201 

M  =  N 

DO  67  K=1,N 
AP(K )=A(M) 
67  M=M-1 

.  DO  432  K=l,30 
ROOTI (K)=0. 
432  ROOTR(K)=0. 
NX  =  NO 

DO  15  <=1,31 
IP  *{Kj=8,fi 
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CALL  POLYRT  ( AP , X , NX , ROOTR , ROOT  I , 1 . E-5 ) 
PRINT  69 

PRINT  68» (ROOTR(K) ,K=1 *NO) 
PRINT  70 

PRINT  68, (ROOTI (K) ,K=1 »NO) 

CALL  DRAW ( NO » ROOTR* ROOT  I »MOD»l »LAB, I T I TLE»XSCALE» YSCALE  ♦ 
11. 7, 2, 2*8,9,0, LAST) 
MOD  =  2 
202  FORMAT  (//////, 21H  THE  SYSTEM  ZEROS  ARE,/) 
K  =  l 

3  IF(B(K)  )  1,2,1 
2  K  =  K  +  1 

GO  TO  3 
1  NORD=N-K 

IF(NORD-l)  6,4,5 

4  ZFRO=-B(K+l )/B(K) 

7  FORMAT  (/,16HTHE  SYSTEM  ZERO= , F15 . 5 , / ) 
PRINT  7, ZERO 

ROOTM( 1 )=ZERO 

ROOTM( 2)=16.*XSCALE 

ROOTJ(1)=0. 

ROOTJ( 2)=0. 

NORD=2 

GO  TO  11 
6  PRINT  9 
9  FORMAT  (/,25HALL  ZEROS  ARE  AT  INFINITY) 

GO  TO  8 

5  NN=NORD+l 

DO  10  L=1,NN 
R(L)=B(K) 

10  K=K+1 
PRINT  202 
M  =  NN 

DO  46  K=1,NN 
AZ (K)=R(M) 
46  M=M-1 

DO  433  K=l,30 
ROOTM(K)=0. 
433  ROOTJ(K)=0. 
NORX=NORD 
DO  16  K=l,31 
16  X(K)=0.0 

CALL  POLYRT ( AZ , X , NORX , ROOTM , ROOT J , 1 . E-5 ) 

PRINT  69 

PRINT  68,  (ROOTM(K)  ,K=1  ,NORD) 

PRINT  70 

PRINT  68,  (ROOTJ(K) ,K  =  1 ,NORD) 

11  CALL  DRAW(NORD, ROOTM, ROOTJ, MOD, 3, LAB, I T I TLE , XSCALE , YSCALE* 
11, 7, 2, 2, 8, 9,0, LAST) 

8  CONTINUE 

GO  TO (3 1,3 2, 3 3, 34, 3 5, 36, 37, 38, 39, 40 ) ,ND 

31  G=1.C076 
GO  TO  41 

32  G=1.016 
GO  TO  41 
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GO  TO  41 
34  6«1 ,031 ? 

00  TO  41 
3C,  0  =  1.0394 

GO  TO  41 
3  6  0=1.0483 

GO  TO  4] 

37  G=1.0568 
GO  TO  41 

38  G=1.0633 
GO  TO  41 

39  G=1.071 
GO  TO  4  1 

40  0=1 .078 

41  CONTINUE 

00  101  J=l,30' 

DO  100  K=l,10 

DO  300  L=ltN 
300  R(L )=A(L)+B(L )*VAR 

CALL  ROOTX 

DO  71  JJ=1,N0 

IF( AB5F (ROOT  I <JJ))-5.E-04)  61  »61  »62 
61  ROOTJ(JJ)=0. 

GO  TO  71 
6?  ROOTJf  JJ)  =ROOTI  (  J  J  j 
71  CONTINUE 

100  VAR=VAP*G 
IF(J-30)  ini,98f98 

98  M0D=3 

101  CALL  DRAW(N0»R00TR,R00TJ,M0D»2»LAB» I T I T LE ♦ XSC ALE ♦ YSC ALE  ♦ 
11»7.2.2»8»9,0.LAST ) 

GO  TO  206 
END 

SUBROUTINE  ROOTX 

DIMENSION  C(31)»D(29)tR(129) »ROOTR( 128) .ROOT I(128).EE(31).FE<31) 
1  ,AP(  129 ) ,X(  129) 
COMMON  R,VAR,NO.ROOTR»ROOT I »AP,X 
M=l 

pO  1492  MA=1»31 

EF(MA)=n. 

1492  FF(MA)=0. 

20  BFTAN=ROOTI (M)+EE(M) 

ALFAN=ROOTR(M)+FE(M) 

DO  7  1=1  »100 

S  =  2.*/\LFAN 

T=-( ALFAN**2+BETAN**2) 

C(  1  )=R(  1  ) 

C(2 )=R( 2)+S*R(  1  ) 

NC=N0+1 

DO  2  L=^»NC 
2  C(L)  =  R(L)+S*C(L-1  )  +  T*C(L"-2) 

AN=  C(N0+1)-ALFAN*C(N0) 

BN=  BETAN*C(NO) 

IF     (NO-3)     21»     lit     18 
17    CN    =    3.*R(  1)*( ALFAN**2-BETAN**2)    +    2 .*R ( 2 ) *ALFAfi  +    R(3) 

$N    =    G.*f3  (  ;  )*A|_FAN*BFTAN    *    2  •  *B  t  2  )  *§£T*N 
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GO  TO  19 
21  CN  =  2.*R( 1 )*ALFAN  +  R(2) 
DN   =  2.*R( 1 )*BETAN 
GO  TO  19 

18  D(l )  =  C(  1) 
D(2)=C( 2)+S*D( 1  ) 
NU=NO-l 

DO  3  N=3.NU 

3  D ( N ) =C ( N ) +S*D ( N-l ) +  T*D ( N-2  ) 
CN=  C(NO)-2.*D(NO-2)*BETAN**2 

DN=  2.*BETAN*(D(NO-l ) -ALFAN*D ( NO-2 ) ) 

19  ALFA=ALFAN-(AN*CN+BN*DN)/(CN**2+DN**2 ) 
BETA  =BETAN  +  ( AN*DN-BN*CN ) / ( CN**2+DN**2 ) 
AZETA=.01+ABSF(BETAN) 

ZETA=SORTF( .01-T) 

IF(ABSF( (ALFA-ALFAN)/ZETA     )-. 1 ) 55 » 55 , 52 
5  5  IF(ABSF( (BET A-BETAN)/AZETA)-1. ) 56,56,52 

56    EF(M)     =(ALFA-ALFAN)*10.    +    BETA    -    ROOTKM) 

FF(M)=(BETA    -    BETAN)*10.    +ALFA    -    ROOTR(M) 

IF(A3SF( ( ALFA-ALFAN) /(ALFAN+1. ) ) -5 . E-4 ) 4 ,4 , 5 

4  IF(ABSF( (BETA-BETAN) /(BETAN+1. ) ) -5 .E-4 ) 6 ,6 » 5 

5  ALFAN=ALFA 
7  BETAN=BETA 

52  CONTINUE 

MR  =  NC 

DO  6777  K=1»NC 

AP(K)=R(MR) 
6777  MR=MR-1 

DO  434  K=1.30 

ROOTI (K)=0. 
434  ROOTR(K)=0. 

NX  =  NO 

DO  83  K=1.31 
83  X(K)=0.0 

CALL  POLYRT  ( AP , X , NX , ROOTR , ROOT  I , 1 . E-5 ) 

GO  TO  16 

6  ROOTR(M)=ALFA 
1  ROOTI (M)=BETA 

12  IF(M-NO)  15,16,16 

15  M=M+1 

GO  TO  20 

16  RETURN 
END 
END 
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CHARACTERISTIC  EQUATIONS 
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